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We discuss the connection between the Kolmogorov-Sinai entropy, has, and the production rate 
of the coarse grained Gibbs entropy, re- Detailed numerical computations show that the (often 
accepted) identification of the two quantities does not hold in systems with intermittent behavior 
and/or very different characteristic times and in systems presenting pseudo-chaos. The basic reason 
of this fact is in the asymptotic (with respect to time) nature of Iiks, while re is a quantity related 
to short time features of a system. 
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I. INTRODUCTION 

The physical entropy is a quantity that plays a key 
role in the understanding of the basic laws ruling the 
macroscopic behavior of systems with many degrees of 
freedom. We just mention the Boltzmann's microscopic 
interpretation of the macroscopic Clausius equilibrium 
entropy, and the celebrated H-theorem 0. 

On the other hand the term entropy is widely used also 
in contexts different from thermodynamics and statistical 
mechanics. In the information theory there is the Shan- 
non entropy 0- 0- 01- while in dynamical systems one 
uses the Kolmogorov-Sinai entropy (and other entropic 
quantities like, e.g., the Renyi entropies) 0,IH- 

The connection between some properties of the non- 
equilibrium thermodynamical systems and the under- 
lying chaotic dynamics, recently, has attracted the in- 
terest of many scientists H, 0,H 0- The main ques- 
tion is if and how the dynamical characteristic quanti- 
ties (such as Lyapunov exponents or Kolmogorov-Sinai 
entropy) are related to macroscopic physical properties 
(e.g.: diffusion coefficients and entropy production rate) 

liiliiiiimmnaninip. 

Some authors, on the basis of reasonable arguments 
and numerical computations on simple dynamical sys- 
tems, claim the existence of a relation between the 
Kolmogorov-Sinai entropy, hxs, and the production rate, 
re of a suitably averaged coarse-grained Gibbs entropy 



The main goal of this work is to show that this is not 
the generic case. To this end we will study two different 
kinds of discrete time systems: 

a) one-dimensional intermittent maps, with local Lya- 
punov exponents very different from the mean 
value; 

b) slightly coupled maps with very different (uncou- 



pled) Lyapunov exponents. 



We will show that when there are different (local Lya- 
punov exponents) time scales in the different regions of 
the phase space, there may be no room for an identifi- 
cation, on a meaningful time interval, between hxs an d 
tq. This is so because, when averaging at a fixed time 
the entropy contributions originating from regions with 
different time-scales, densities with different space-scales 
get mixed, at possibly different relaxation stages. 

A second point we want to stress is that, even if a 
linearly increasing time behavior of the coarse-grained 
entropy is observed, the rate of growth is not necessarily 
given by the Kolmogorov-Sinai invariant of the system. 
This will be shown by studying discretized versions of a 
chaotic map (i.e. an automaton). In such a case, the 
memory of the chaotic character of the original system 
may allow for a transient chaotic-like regime long enough 
to make the behavior of the strictly periodic system in- 
distinguishable from the one of its truly chaotic ancestor. 
That is, pseudochaos is at work 0,0]: the long time 
properties of the system, like hxs{— 0), remain hidden 
and do not affect quantities, like tq{^ 0), related to short 
time features of the system. 

The paper is organized as follows. In Sect. 2 we recall 
some basic concepts and methods in chaotic dynamical 
systems and statistical mechanics. In addition we give a 
simple argument for the connection between hxs and tq, 
stressing the weak points of the argument. In Sect. 3 we 
discuss the results of numerical computation for the time 
evolution of the coarse-grained Gibbs entropy in systems 
with "non trivial" dynamical features, i.e. with intcr- 
mittency and different characteristic times. Sect. 4 is 
devoted to conclusions and discussions. 
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II. A BRIEF OVERVIEW OF BASIC FACTS 

Since the pioneering work of Kolmogorov 0, the 
relevance of the Kolmogorov-Sinai entropy for a proper 
characterization of the behavior of a chaotic dynamical 
system was clear. To perform the computation of Kks 
one has to choose a partition A of the phase space and 
to assign each cell of the partition an integer value i. In 
such a way, the trajectories of a dynamical system, with 
continuous states, sampled at discrete times, 

x(l),x(2),...,x(j),...,x(T) 

become symbolic sequences 

i(l),i(2),...,i(j),...,i(T) (1) 

whose meaning is that at time j the trajectory is 
in the cell labeled by Then one defines the 

probability of each word (or block) of length n, 
p^'(ki, k%, .., k n ), counting how many times one meets 
the word k\, hi, fc„ along the sequence {Q. The en- 
tropy of the blocks of size n, H^'(n), thus reads 

H {A \n) = - Y, P {A) {k u ..,k n )\ogp^\k u ..,k n ). 

ki , . .,k n 

(2) 

The Kolmogorov-Sinai entropy is defined by the sup over 
all partitions of the asymptotic value of the rate of in- 
crease of H^(n), i.e. 

hxs — sup lim . (3) 

_4 n^oo n 

The quantity hxs is a numerical invariant that gives a 
good characterization of a chaotic system, but unfortu- 
nately it is almost impossible to compute analytically 
(except very few simple cases) and also rather difficult 
from a numerical point a view. From a physical point 
of view it is rather natural to use regular partitions with 
hyper-cubic cells of edge e. Let us denote with H^ e \n) 
the n-block entropy on a partition of this kind and e- 
entropy the limit 

h( e ) = lim g(c) (") = lim H (t \n+l)-H (e) {n) . (4) 

n — >oo 77, n — >oo 

It is a remarkable fact that h(e), computed with different 
values of e, can give very interesting information about 
the properties of the system^, 0. Moreover it is pos- 
sible to obtain h^s by considering the limit e — > in 
(0J , instead of the sup operation in J2J • By recalling the 
Shannon-McMillan equipartition theorem (stating that, 
for large n, the number of "typical" n-words increases 
as exp[/i(e) n]) one has that h{e) gives the (asymptotic in 
time) exponentially growing rate of the number of typical 
trajectories of the system, in the limit of high resolution 
(as measured by e). It is intuitive that this growing rate 
must be linked to the exponentially fast separation of 



nearby trajectories. The Pesin theorem is the rigorous 
statement of this idea: all the expanding directions con- 
tribute to the diversification of the trajectories and to 
the increase of their number 

Ai>0 

where (A^ > 0) means the sum over positive A^. Pcsin's 
identity JSJ provides us with a useful alternative way 
to compute \iks- The Lyapunov exponents can be nu- 
merically computed without particular difficulties, even 
in high-dimensional systems; on the contrary, because 
of the exponential proliferation of the n-words, the 
Kolmogorov-Sinai entropy becomes rapidly unattainable 
by numerical methods (a part low-dimensional systems). 
Therefore very often the Pesin formula is basically the 
unique way to compute hxs- 

Notice that hxs is an entropy rate defined on the en- 
semble of the trajectories of a system (according to some 
stationary probability measure), even if, by means of er- 
godicity, it is practically computed using only one single 
long trajectory. Qualitative arguments may be given [ljj 
to support a connection between hxs and the rate of 
variation of an entropy-like quantity, that we call Gibbs 
entropy, defined on the phase space of the system. A 
possible line of reasoning is the following. 

Consider a deterministic dynamical law 

x -> T*x (6) 

(where x is a £>-dimensional vector) and a probability 
density p(x,t), that gives a distribution of states of the 
system throughout its phase space at a time t. We define 
the Gibbs entropy of p as follows 

S( Pt ) = - J p(x, t) ln[p(x, t)]dx , (7) 

i.e. its conditional entropy with respect to a uniform 
density. For chaotic dissipative systems, where p(x, t) 
tends to a singular (fractal) measure, definition @ be- 
comes meaningless. Nevertheless, following a nice idea 
of Ruelle ^2$, one can avoid this difficulty, simply by 
adding (or considering unavoidably present) a small noise 
term in the evolution law: in such a way one obtains a 
/o(x, t) continuous with respect to the Lebesgue measure. 
If J(x.,t) is the Jacobian of ©, a straightforward com- 
putation gives 

S(pt)=S(p ) + J p(x,t)ki\J(x,t)\dx . (8) 

In the case of volume conserving evolutions, one has: 
S(pt) = S(po). To allow for an entropy variation one 
needs a coarse-graining. Let us consider a hyper-cubic 
partition, as introduced above, and let us define the prob- 
ability P e (i,t) to find the state of the system in the cell 
i at time t: 

P t {i,t)= f p(x,t)dx (9) 

J At 
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where Af is the region singled out by the i—th cell. Let 
us introduce the e-coarse-graincd Gibbs entropy 



If e is small enough S(pt) and S e (Pt) are trivially related: 

r 



S e (P t ) ~ S(p t ) + Dhx 



If one considers a distribution of initial conditions that 
is different from zero only over one (or very few) cell(s), 
one has, for e small enough and a time not too short, 



5 e (P t ) = S e (P ) + h KS t. 



(10) 



One can argue as follows. Assume that the system has 
m positive Lyapunov exponents and that p(x, 0) is local- 
ized around x c (0). In a suitable reference system (with 
the axes along the eigendirections of the Lyapunov expo- 
nents), if p(x, 0) has a Gaussian shape, p(x.,t), for some 
times, is still well approximated by a Gaussian with vari- 
ances 



o-|(i) = cr] (0) exp{2A^} 



therefore: 



D 



J= i */27raf (t) 



(11) 



(12) 



where x c (i) is the state evolved from x c (0). From this, 
in the non-grained case, one gets 

s( Pt ) = s(p ) + E ln ^# = s m + E A ^ • 

It is clear that S(p t ) = S(po) if the phase space volume is 
conserved. Considering now the coarse-graining @, one 
has that along the directions of the negative Lyapunov 
exponents (m + 1, m + 2, . . .), for a long enough t: 

o*(t) ~ o-fc(0)e-' A *l* < e 
This implies that 



™ * n 



i 



i= i yj2na?(t) 



and therefore 



5 e (P f ) = ^(^o) + EV- 

With the aid of the Pesin's formula eq. I|1U|) follows. 
Let us stress that the transition from JSJ to (|10|) is allowed 
by the fact that, in the presence of a coarse-graining, the 



contracting eigendirections (corresponding to the nega- 
tive values of the Lyapunov exponents) cannot balance 
the effects of the expanding ones. 

At this point we have to notice that, by definition, the 
Gibbs entropy J7J explicitly depends on the particular 
chosen initial density. In the discussion here above this 
dependence may be labeled by the cell, x c , where the dis- 
tribution is initially different from zero. On the contrary, 
fix s is an asymptotic global property of the system. So, 
one may expect that a density independent behavior, as 
in (|10f> . can be found only in "friendly" dynamical sys- 
tems, i.e. systems with no fluctuations. A generic system 
possesses a certain degree of intermittency, so that, for 
instance, the expanding and contracting properties may 
strongly depend on the phase space region the trajectory 
is visiting. This calls for an averaging over the initial con- 
dition x c , weighted, say, by the natural invariant measure 
of the system: 



S(pt) - S(t) 



J S(t\x c 



lp e ,(x C )dx C 



(13) 



where S(t\xc) is S(pt) with po( x ) localized around x c . 
The same averaging procedure leads to S e (t) from S e (Pt). 
This operation yields intrinsic quantities, that can de- 
pend on global properties of the system. The interesting 
question is whether the simple relation (|10|) survives, as 
an observable property, for the averaged coarse-grained 
entropy Notice that equation (JTDJ , besides its concep- 
tual interest, would result in a numerically simple way to 
determine h^s- I n some systems suitable ranges of t and 
e exist where the relation is verified |16j . However, we 
believe, in agreement with Ref. that what has been 
found in Ref. 0] is just a "lucky" coincidence. With 
regard to that, it is important to stress that in the above 
presented arguments, for the derivation of l|l(JII . there are 
(at least) two delicate points: 

a) both Lyapunov exponents and Kolmogorov-Sinai 
entropy are quantities defined in the limits of high 
resolution (e — > 0) and long times (t — > oo); 

b) a behavior like (|11|) holds only for "short" time, 
i.e.: 

^aT 111 ^) ; 

and this is so also for the linear behavior in 11U|) 
that is expected to be valid for 

t<W- 

Ai e 

If intermittency is present one has to replace Ai 
with the largest local Lyapunov exponent. 

Since for a coarse-grained Gibbs entropy a finite e is 
mandatory, it is not obvious that the previous time 
regimes a) and b) have a non-empty overlap. One may 
note that in the entropic analysis of the n-words one of 
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the two asymptotic limits can be relaxed: i.e. one can 
work with non-infinitesimal e and therefore obtain the 
e-entropy h(e) which is an asymptotic in time quantity 
associated with a finite tolerance e. However, once the 
size of the cells for the coarse-grained Gibbs entropy is 
fixed, nonetheless p(x, t) evolves developing structures on 
scales l(t) w eexp(At) increasing in time. Therefore it is 
not trivial at all that a simple relation between dS e (t)/dt 
and h(e) exists. 

III. A NUMERICAL STUDY OF SIMPLE 
DYNAMICAL SYSTEMS. 

In this section we present an entropic analysis of sim- 
ple dynamical systems which show, in spite of their low 
dimensionality, non trivial features. 

From a practical point of view the computation of S £ (t) 
is performed as follows: 

1. we select several starting conditions xj(0), 34(0), 

x] v (0) all located in the j-th box of size e D 

2. we let evolve all the N 3> 1 starting conditions up 
to a time t obtaining x , 1 (t),x 1 2 (t), ...,x 3 N (t) 

3. we calculate 

1 N 

with 

5(x k) = { 1 if X 6 A * 
v ' ' ) otherwise 

4. we compute the entropy of p e ' J (k,t) defined as 

N 

S V,t) = -$> £J (M)logp £J (M) (14) 
fc=i 

5. we average this quantity on the coarse-grained in- 
variant distribution p e q(j) obtained from /Q e g( x ) 
with the procedure of Eq.@. Thus we have 

N 

S*(t) = J2Pe q (j)S € (j,t). (15) 
J'=l 

A. Intermittent map 

A recent pape r Il5l ] shows that the value of tq in the 
Manneville map [21| significantly differs from the value of 
h ks- These authors perform the analytical calculation of 
S(t) starting from only one particular condition, namely 
the box containing the point x = 0, that in the Man- 
neville map is the point with the lowest local Lyapunov 



exponent. Since in this map the value of the local Lya- 
punov exponent ranges from very low values (A(a;) — > 
for x — > 0) to values considerably grater than 1, the au- 
thors correctly argue that the discrepancy is due to the 
variability of \{x). Nevertheless the authors of do 
not perform the final average over the initial condition. 
In the following we see that the averaging procedure i|15|) 
is not able to recover the condition of Eq. I|1(J|I neither in 
the Manneville map nor in a less-intermittent map. Let 
us start with a simple case, i.e. the modified tent map 
given by 

_ / X t/P if x t <p 

Xt+1 ((1- xt)/(l -p) if p < x t < 1. [Lb) 

This map is the standard tent map if p — 1/2, while for 
small values of p one has an intermittent behavior charac- 
terized by two very different local Lyapunov exponents, 
namely: 

A+ = -logp, A_ =-log(l-p). (17) 

The stationary distribution is constant between and 1 
and consequently we have 

Iiks = pX++(l-p)X- = -plogp-(l-p)log(l-p) (18) 

assuming its maximum value for p = 1/2 and decreasing 
its value for p — > 0. Note that this map presents an 
intermittent behavior, but not in a critical way like in the 
Manneville's map, and the distribution of the length of 
the "laminar" zones (the permanence in the zone [p, 1]) is 
simply exponential without power law tail. Nevertheless 
we see that already in this system Eq. (|l(J|) does not hold. 

The numerical results are shown in fig^ We set 
e = 10~ 3 and we scale the time axes with the inverse of 
hxsip), given by Eq. l(T%|) . so if Ea. l(TU|l holds we should 
observe that all the curves, with different p, collapse on 
the straight line S e (t) = fiKs(p)t- As one can see the 
agreement is good only for p = 1/2 (with no intermit- 
tency) while becomes worse and worse decreasing p. The 
main point is that the linear behavior of S e (t) should 
hold, according to the heuristic arguments of Ref. [l6| . 
till to a time given by 

t Hn ~ log (^j (19) 

corresponding, on the scaled time of Fig^ to the value 
— loge ~ 7. This is observed for p = 1/2 but for small 
values of p one has that S e (t) increases in time with a 
"wrong" slope (different from hxs) an d later it exhibits 
a rather long crossover. 

The origin of this effect is in the intermittent behavior 
of the system. Indeed the realizations S e (j,t) starting in 
the zone [Q,p] are spread on the whole interval [0, 1] after 
few steps almost reaching the asymptotic value of — log e 
while the "realizations" starting for example near the un- 
stable equilibrium point x = 1/(2 — p) takes several time 
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steps to reach the saturation giving a dominant contribu- 
tion to the rate of increase of S e (t). In this way the S e (t) 
computed with Ea. H15(l does not increase in time follow- 
ing the naive argument yielding to Eq. l|l(Jf> . The reason 
of the discontinuity in S f ~(l) for p — > 1/2 is explained in 
Appendix A. 

Fig[21shows other interesting properties of the behavior 
of S e (t) for p — 0.1 and p = 0.4 at varying the value 
of e. As one can see, in the slightly intermittent case 
p = 0.4 the rescaled curves collapse together confirming 
the assumption of Ea. (|10f> : while in the intermittent case 
p = 0.1 the curves do not collapse and only for very low 
values of e a linear growth of 5 e (t) is present. Let us 
stress the fact that, at variance with the case shown in 
Fig. [2)3, in the intermittent case (FigOJO the crossover 
regime (after the linear one and before the saturation) 
is very long and is comparable with the duration of the 
linear regime. 

To connect our results to the ones of Ref. [H wc also 
perform the same calculation for the Manneville map 
given by 



x t +i 



x t + kx* 
(1 -**)/(! 



d) 



x t < d 
d < x t < 



1. 



(20) 



where d fulfills 



d + kd z = 1. 

In the range 3/2 < z < 2 and k <C 1 there is an invariant 
distribution [4J, Il5| . Under these conditions the distribu- 
tion of the permanence time r has a power law tail with 
(r) < oo while (r 2 ) diverges. The only difference from 
the tent map calculation is that the invariant distribution 
p eq (x) is not constant and therefore we have to compute 
it numerically. Also in this case we observe a "wrong" 
slope (lower than hxs) and an extremely long crossover 
behavior. 



CO 




p = 0.500 
p = 0.240 
p = 0.120 
p = 0.060 
p = 0.030 
p = 0.015 
t h K s(P) 



3 4 5 6 
t h KS (P) 



10 



FIG. 1: S e {t) as a function oith K s(p) with e = 10" 3 and dif- 
ferent values of p. Note that Eq. llOt should give the straight 
line thics(p) for all values of p. 
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FIG. 2: S E (t)/S c (oo) as a function of t/S E (oo) with different 
values of e. In fig. a) p = 0.1, in fig. b) p = 0.4 



B. 2D map 

Let us now discuss a system with two different time- 
scales. We choose the following 2D maps for the variable 
x t and y t 



x t +i= f(x t ) + crcos(27r(x t + y t )) 
yt+i = ry t + cr cos(27r(x t + y±)) mod 1 



(21) 



with 



f(x) 



if x < p 
if p < x < 1 



(22) 



x/p 

(l-x)/(l-p) 

and r is an integer greater than 1. The map for x t is a 
modified tent map, like in the previous section, while the 
map for yt is a generalized Bernoulli shift. We choose 
p = 0.3 and r = 5 in order to have two largely different 
uncoupled Lyapunov exponents, namely ~ 0.61 and 
X y = log 5 ~ 1.61. We use this kind of coupling to avoid 
discontinuity at x — or x = 1. 

We will see that in this system the two Lyapunov ex- 
ponents X x and X y have a role rather similar to that of 
A + and A_ for the system (fTp|l . 

The coarse-grained Gibbs entropy of this 2D map is 
calculated starting from a 2-dimensional cell of linear 
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size e and performing the average over the two dimen- 
sional coarse-grained invariant distribution, that also in 
this case turns out to be flat in the slightly coupled case. 
We also numerically calculated the Lyapunov exponents 
in the coupled case: for the small values of a used in the 
numerical calculations, we observe no relevant changes 
from the uncoupled case. We set e = 5 ■ 10~ 3 and we 
study the behavior of S e (t) with different values of a, 
in such a way we can study the crossover behavior at 
varying a. The results for are shown in fig El 
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FIG. 3: S*(t) as a function of t with different values of a, 
e = 5- 10" 3 , p= 0.3, r = 5 

As one can see, at very small t, S e (t) increases with 
a slope equal to the Kolmogorov-Sinai entropy hxs — 
X x + X y = —plogp — (1 — p) log(l — p) + logr. After that 
the phase space of the variable y saturates and the slope 
changes. In the case a < e the influence of the coupling 
is almost negligible and the second slope is the same we 
should observe in the ID case. In other words, as shown 
in the figure, we have S e (t) ~ S^(oo) + S x (t) where S x (t) 
is the entropy calculated from the marginal probability 
density p x {x, t) — J dyp(x, y, t). In a similar way Sy(t) is 
obtained from p y {y, t) = J dxp(x, y, t). 

When a w e the noise on the slowest variable x is 
enough to spread the distribution as fast as it happens 
for the variable y so no change in slope is observed. 

Thus we observe that, for values of a > e, tq cor- 
responds quite well to the value of has in a broad 
time interval while for values a <^ e this correspon- 
dence longs only for a very short time. This is note- 
worthy because this map has pratically the same Iiks 
entropy at the different values of a used (for a = 
h KS = -plog(p) - (I - p)log(I -p) + log(r) = 2.22 
while for a — 0.1 we numerically obtain fixs — 2.19) but 
with values of a > e we have tq ks h^s while for a < e 
the equivalence is lost or longs for a very short time. 

It is interesting to compare Fig|3|with FigJ2Ju: in both 
cases the "naive" behavior (i.e. S e (t) — S e (0) = hxst) 
may have a very short duration and there is a long 
crossover. The origin of this crossover is due to a sort of 
"contamination effect" of different times (in other words 



different mechanisms) involved. For the system (116(1 this 
is due to intermittency (the system "feels" A+ and A_) 
while for the 2D map it is due to the existence of differ- 
ent times that are relevant at different spatial resolution 
scales. 

In systems with many degrees of freedom the "contam- 
ination effect" can produce rather impressive behaviors. 
As example we can cite the case of fully developed tur- 
bulence where, because of the existence of many different 
characteristic times, the growth of the distance between 
two trajectories is a power law in time instead to be an 
exponential (although the system is chaotic) |2^| . 



C. Discretized ID map 

We now discuss a class of dynamical system with zero 
Hks'- the discretized maps also called coarse-grained de- 
terministic automata j24[. We study the discretized (in 
phase space) version of ID map defined by the equation 



n t +i = ?/ 



/(«*) 



(23) 



where r\ is the discretization parameter (the number of 
discretized state is l/?7), [-J denotes the integer part and 
fix) is a tent map: 

{x I p if x *C p 

(l-x)/(l-p) if P <x<l. ( 24 ) 

Even if this system is the discretized version of a 
chaotic map the finiteness of the available states forces 
the dynamics to become periodic. In ref.|24j it has been 
observed that the block entropy H e (n) defined as in 
Eq.(0| with e ^> r\ (i.e. several discrete states in each 
cell) for small values of n increases with a slope given by 
the Kolmogorov entropy Hks of the corresponding con- 
tinuous system Xt+\ — f(xt) while for 



>t, 



1 



h K S 



log?7 



(25) 



the block entropy stops increasing thus revealing the pe- 
riodic nature of the dynamics. We present the behavior 
of S e (t) in this system, varying the value of p in order to 
observe also the effect of intermittency. As figQ] shows, 
the coarse-grained Gibbs'entropy increases as a function 
of the scaled time tfiKs{p) with a slope comparable with 
1; while it reaches the saturation values at scaled time 
given approximatively by 



t. 



-r— log ( - 

h K s \ e 



7. 



Note that, in spite of the fact that the KS-entropy of 
(|23|l is strictly zero, this behavior is practically the same 
of the continous system, thus revealing one of the main 
problem of the use of the coarse-grained Gibbs'entropy 
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to detect Hks- Indeed the block entropy H e (n) "feels" 
that the system is periodic only for n > t p , but 

hicstp — ~ logr? > - log e. (26) 

The last inequality follows from the physical condition we 
use i.e. that e ^> rj in order to have several state inside 
each cell. The main point is that since e 3> i] 

t p - t — log(l/r?) 
h K s 

is larger than the saturation time t sat . Therefore the 
coarse-grained Gibbs 'entropy saturates well before the 
system feels to be periodic. 
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FIG. 4: S e (t) as a function of thKs{p) with different values 
of Pl e= 10" 3 , 1/t] = 3- 10 6 



IV. CONCLUSIONS 

We have shown that, in spite of some folklore, there 
is a rather loose relation between the Kolmogorov-Sinai 
entropy and the growth of the coarse grained (Gibbs like) 
entropy Such claimed connection exists only in very spe- 
cial cases, namely systems with a unique characteristic 
time and very weak intermittency (i.e. small fluctau- 
tions of the local Lyapunov exponent). On the contrary, 
in more interesting (and closer to reality) systems, with 
multiple characteristic times and/or non negligible inter- 
mittency, the relation between hxs an d tq holds (if any) 
only for a very short time. 

The main reason of this is due to the asymptotic nature 
of hfcs (as well as the e-entropy), i.e. its relevance at 
very large time intervals. On the contrary the growth 
of the coarse grained entropy only involves short time 
intervals, and during the early time evolution of 5* e (t) one 
has entanglement of behaviours at different characteristic 
space scales. 

This phenomenon is rather similar to that observed in 
the spreading of passive tracers in closed basins [25|, In 



such a case, if the characteristic length scale of the Eule- 
rian velocities is not very small, compared with the size 
of the basin, both the diffusion coefficient and the Lya- 
punov exponent do not give relevant information about 
the mechanism of spreading. 

We stress again that the failure of the relation between 
Hks and tq is due to the fact that the growth of the 
coarse grained entropy is not related to asymptotic prop- 
erties (i.e. large time and small resolution). This is par- 
ticularly evident in deterministic discrete states systems 
that, although non chaotic (i.e. with hxs = 0), show a 
behaviour of S e (t) very similar to that observed in gen- 
uine chaotic systems (with hxs 7^ 0). 
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APPENDIX: A 

We discuss here the discontinuity in S e (l) for p ^ 
1/2 observed in fig^ Indeed for p = 1/2 we have 
5 e (l) = log 2 while for p = 1/2 - £, with £ <C 1, we 
have 5 e (l) ~ 0.85. Obviously S" e (0) = for all values of 
p. To understand the reason of this discontinuity let us 
compute the value of 5 e (l). If we start from the j-th cell 
we have 

P ed (k,0) = S kJ . 

If p = 1/2 at time 1 this cell will be spread on exactly 
two cells and the trajectory starting at the border of a 
cell will go exactly on the border of another cell so if, for 
sake of simplicity, j < l/2e (i.e. if the cell is in the region 
x < 1/2), we obtain 

P e,3 (k,l) = ^(S k ,2j + S k ,2j+l) 

leading to S e ( 1 ) = log 2 for each j . If we have p = 1 / 2 — £ 
with ( < 1 the size of the phase space region populated 
after one time step is essentially the same. Neverthe- 
less this time the trajectory starting at the border of a 
cell does not always go to another cell border so in the 
averaging procedure of Eq. (|15fl there will be some start- 
ing conditions giving approximately a superposition of 
the cell borders after one time step (for example i = 0) 
and consequently a value of S' £ (l) log 2, but also some 
starting conditions j resulting in 

p e ' 3 {k,l) = ~4,2j+l + ^(4,2j +5 k ,2j+2). 

These terms give a contribution to the average greater 
than log 2 and lead to the discontinuity of fig^] when 
going from p = 1/2 to p ^ 1/2. Anyway this behavior 
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is not important for greater values of time so it does not 
influence the behavior of S e (t) for t > 1. 

In order to check our numerical results we also com- 
pute S e (t) from p(x, t) obtained by the Perron- Frobenius 
equation. In the case of the tent map we use directly the 
expression of the formal evolution of p(x, t) given by 

p{x, t + 1) = pp{px, t) + (1 - P )p(l - x(l -p),t) (A.l) 

implemented in a recursive algorithm on a computer pro- 



gram. We fix 

p(xO) = i 1/€ ifxeA ^ 
Hy ' ; [0 otherwise 

we let p(x,t) evolves and then we compute S e (t). The 
results obtained with this method are in perfect agree- 
ment with those in Section III and confirm the presence 
of the discontinuity. 
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